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1. Introduction 

There exist many ways one can choose to approach the subject of these lectures. Taken that 
the presentation is hmited in time (spent by the students and lecturers at the Caribbean seashore 
during the Mochima school in theoretical physics) and in space (given by the Editors of this 
volume) this choice becomes a difficult one. The way of presentation I decided to follow was chosen 
for several reasons. Extremely high theoretical level of certain modern renormalization group (RG) 
studies of criticality in different systems (i.e. the "language" of these studies) sometimes does not 
allow uninitiated reader to follow the derivations and even to understand the problem statement 
and/or physical consequences of the results. Therefore the goal would be to provide a minimal 
vocabulary, explaining main notions as simple as possible. On the other hand, it is tempting to 
use such a simple vocabulary to make a short review of state-of-the art RG studies in a certain 
domain. RG explanation of criticality in "non-ideal" 3d systems might be a good candidate for such 
a domain: it is a subject of ongoing activity where important results have recently been obtained 
and still a lot is to be done. 

The following account will serve this purpose: after mentioning several examples of criticality 
and scaling in condensed matter physics (section [5J) and introducing model Hamiltonians of the 
"non- ideal" systems we shall be interested in (section O I shall give the main ideas and notions 
of the renormalization taking as an example a simple Id Ising model (section 0)). Once the reader 
is acquainted with the RG transformation, its flow and fixed points, stability, universality and 
scaling, I shall pass to the "non-ideal" systems showing how to obtain their effective Hamiltonians 
and to the reviewing of recent results in this domain (sections and EJ. Some conclusions and 
outlook are given in section [7| 

2. Criticality and scaling 

It is generally recognized that the term "critical point" was introduced in 1869 by Thomas An- 
drews who studied a special point for carbon dioxide at about 31°C and 73 atmospheres pressure 
where the properties of liquid and of gas become indistinguishable [ 01 • Approaching the critical 
point Tc from below, the liquid-gas density difference obeys a power law scaling governed by the 
critical exponent [3: 

PL-pG^{Tc-Tf, T^T-, (2.1) 
as shown in the figure Similar power laws were found for some other physical quantities. 




Figure 1. Order parameters for the critical point of fluids (a) and Curie point of magnets (b). 

describing their approach to zero or singularity at Tc- Less known is the fact, that the mean 
field values of critical exponents (which are integers or their simple relations, e.g. /3 = 1/2 for 
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equation (|2.1|l l were questioned already over a century ago [El^- Therefore much time passed 
before there appeared the theory that both explained scaling behaviour 1)2. l|l found in condensed 
matter physics, and, on a larger scale, not only in physics, as well as offered quantitatively accurate 
predictions for the exponents. As we now know this theory is based on the renormalization group 
formalism and has its origin in quantum field problems 

Next experimental observations of critical behaviour arrived in 1895, when Pierre Curie showed 
that the ferromagnet iron also displayed a special point. This point is the highest temperature at 
which iron can be permanently magnetized at zero external magnetic field which is now called the 
Curie point. Curie himself noticed the parallelism between density-temperature curve at constant 
pressure for carbon dioxide and magnetization-temperature curves at constant magnetic field for 
iron. Behaviour of the order parameter, spontaneous magnetization AI , in the vicinity of Curie 
point (we denote it Tc as well) is governed by the familiar scaling law, cf. equation (|2.1(l : 

MociTc-Tf, T^T-. (2.2) 

It appears that different magnetic systems possess different sets of critical exponents but there 
exist wide classes of magnets characterized by the same values of exponents. For example, uniaxial 
magnets (like FeF2, MnF2, K2C0F4, Rb2CoF4) possess the same critical exponents. Moreover, they 
coincide with those of simple fluids (Xe, SFg, CO2 . . . ) and some other systems. In such a case 
it is said that these systems belong to the same universality class. The reasons for this as well as 
the main characteristics of a system which define the universality classes will be cleared up in the 
subsequent sections. 

As far as the second derivatives of Gibbs free energy (isothermal magnetic susceptibility for 
magnets and isothermal compressibility for fluids) diverge at Tc such phase transitions are referred 
to as the second order phase transitions. This classiflcation is due to Paul Ehrenfest. The order 
parameter does not possess a discontinuity at T^. it is an example of a continuous phase transition. 
Phenomena which occur in the vicinity of the 2nd order phase transition point are called critical 
phenomena. Changes of system structure at Tc are reflected in the behaviour of the pair correla- 
tion function G2{r) (density-density function for fluids and magnetization-magnetization one for 
magnets) and described by two universal exponents 77 and v. The exponential decay of G2{r) with 
distance at T > Tc transforms to the power law one: 

G2(r)^e-'-/«, T>Tc; G2(r) ^ , T = , (2.3) 

where ^ is a correlation length which diverges at Tc'. (, ^ \T ~ Tc\^'^, T — + Tc. 

The XX century witnessed numerous experiments performed after pioneering works of Thomas 
Andrews and Pierre Curie where it was shown that criticality and scaling accompany not only the 
second order thermodynamic phase transitions. They are found in quantum phase transitions, per- 
colation, non-equilibrium dissipative phase transitions. Properties of long flexible polymer chains 
in good solvents are described in terms of critical phenomena as well. This list can be continued. 
However already the above mentioned field of phenomena is very wide, especially if one takes into 
account that any of them occurs in objects differing by their microscopic nature. Thermodynamic 
2nd order phase transitions occur in magnets, ferroelectrics, quantum liquids; percolation phenom- 
ena occurs in resistor networks, random magnets, gas filters. Nevertheless there are at least two 
essential features common for the above mentioned phenomena. These are: the singular character 
of change of properties in certain critical points and universal behaviour in the vicinity of these 
points. As it became clear now the reason for both is anomalous growth of fluctuations and their 
correlation at very high distances in the vicinity of critical points. Correlation range becomes the 
only characteristic scale of the system and this causes the insensitivity of its behaviour to the 
so-called microscopic parameters. In particular this results in the scale invariance of the system 
where critical phenomena occur. 

Since the pioneering work of Ising theoretical description of many-particle systems is often 
based on the so-called classical spin models. They appeared to be of primary importance to reveal 

thank Reinhard Folk for pointing me the work of Jules-Emile Verschaffelt. 
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the main features of criticality. In the next section|31we shall discuss several models currently used 
in describing the critical phenomena, the analysis of which will be a subject of subsequent sections. 

3. Non-ideal spin Hamiltonians: single ion anisotropy, structural disorder, 
frustrations 

In classical spin models, each particle is imitated as a vector located on a lattice site and 
"interacting" with other vectors. The interaction is chosen to mimic the interparticle one whereas 
the dimensionality of the vector is equal to the number of internal degrees of freedom of the particle. 
Below we start with what we call an ideal system, which is a m-vector model and show how to 
include different types of non-idealities within a spin model formalism. 

3.1. An ideal system: m-vector model 

The model describes a system of m-dimensional classical vectors ( "spins" ) located in the sites 
of rf-dimensional hypercubic lattice. The Hamiltonian of the m-vector model in the absence of an 
external magnetic field reads: 

H = ^ J(|R-R'|)5r5r-, (3.1) 

R,R' 

where J(|R — R'|) is an interaction between spins S'r and Sr' located in sites R and R' (we shall 
consider it to be the short-range one and of ferromagnetic origin: J(i?) > 0) and 5r5r' means a 
scalar product. Such a model appears as a natural generalization of the Ising model for the case 
of m-component spin [ 151 and serves as a basic model in describing phase transitions in systems 
with multicomponent order parameter. Below, we shall be mainly interested in the 3d systems. 
However, for the sake of completeness let us recall that ferromagnetic ordering does not occur in 
this model at d = 1 [ B , it does not occur for d = 2 and m > 1 either [ '7 . Whereas the celebrated 
Onsager solution [ 8 of the 2d Ising model (c? = 2, m = 1) brings about the second order phase 
transition with non-trivial values of the critical exponents: (3 — 1/8, v = 1, rj = 1/4, 7 = 7/4, 
a = (exponents 7 and a govern power law scaling of the magnetic susceptibility and specific 
heat, the latter being logarithmically divergent for the 2d Ising model). At d = 3, ferromagnetic 
ordering occurs at any m (see tabled for typical values of exponents), and d = 4 is the upper 
critical dimension of the problem: exponents attain their mean field values for d ^ 4. 

Table 1 . The standard values of the critical exponents of the 3d m-vector model (obtained in [ 
|n| from the high-order renormalization group expansions). 



m 


7 


V 


^ 


/5 


a 


1 


1.2396(13) 


0.6304(13) 


0.0335(250) 


0.3258(14) 


0.109(4) 


2 


1.3169(20) 


0.6703(15) 


0.0354(250) 


0.3470(16) 


-0.011(4) 


3 


1.3895(50) 


0.7073(35) 


0.0355(250) 


0.3662(25) 


-0.122(10) 



We refer to the above model as to the ideal one: it describes the lattice system without any 
defects of structure or complexities of interaction which are often encountered in real systems. Let 
us show how those can be considered within the same formalism. 

3.2. Single-ion anisotropy 

Real substances often are anisotropic. For instance, in cubic crystals one expects the spin 
interaction to react on the lattice structure (crystalline anisotropy) suggesting additional terms in 
the Hamiltonian, invariant under the cubic group [ I10|. Such a single-ion anisotropy breaks the 
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rotational symmetry of the m- vector magnet and the Hamiltonian reads: 

n^^-Y^Jm- R'l) SnSn' + V^E E i^iO' ' (3.2) 

R,R' R i=l 

where V defines anisotropy strength and makes the order parameter to point either along the edges 
{V > 0) or along diagonals {V < 0) of a m-dimensional hypercube. Therefore a model with the 
spin Hamiltonian (|3.2|) is often called a cubic model. 

An interesting phenomenon is observed in 3d cubic magnets: for low values of m their critical 
exponents coincide with those of the "ideal" m-vector model (one speaks about isotropization of 
critical fluctuations), whereas for large m they belong to the new universality class. The value 
m^"*^ at which one regime is changed to the other one is called the marginal dimension. Moreover, 
transition to the low-temperature phase may also occur via the first-order scenario. 

3.3. Structural disorder 

To treat structural disorder within the lattice model one usually introduces random variables 
into the Hamiltonian (|3.1|l and couples them to spin degrees of freedom. The new model can either 
mimic a quenched system (new variables are randomly distributed and fixed in a certain config- 
uration) or an annealed, equilibrium one [lllj. Here, we shall be interested in changes of critical 
behaviour caused by quenched disorder. Furthermore, we shall consider two different examples, 
showing how to introduce disorder via dilution and random anisotropy^. 

3.3.1. Random-site dilution 

To describe dilution, one may introduce the random-site m-vector model, considering situation, 
when spins S'r in occupy only a part of the lattice sites, Ni, the rest N — Ni sites being empty 
(or occupied by non-magnetic atoms). Magnetic and non-magnetic sites are randomly distributed 
and fixed in a certain configuration. The model Hamiltonian reads: 

R,R' 

where cr are the occupation numbers: 

f 1, site R is occupied, 

[ 0, site R IS empty. ^ ' 

Here we do not touch upon the phenomena occuring near the percolation threshold Cpcrc and 
consider the so-called weak dilution, c = Ni/N ^ Cpcrc. To complete the model one should choose 
a certain distribution function for the occupation numbers cr. Let us consider the case when the 
site is occupied with the probability c and is empty with the probability (1 — c) and this probability 
does not depend on the occupation numbers on the neighbouring sites. Such a situation corresponds 
to the following occupation probability V{{c^}): 

^({cr}) = 11^(^1^)' p(cr) = c(5(cr-1) + (1-c),5(cr). (3.5) 

R 

For different m, the above model describes magnetic phase transitions in crystalline alloys of uniax- 
ial magnets and their non-magnetic isomorphs FexZni_xF2, MnxZni_xF2 [Eli diluted Heisenberg- 
like magnets [ 1141 115| , supcrfluid phase transition in He^ in porous medium [ I16j . Similar to the 
cubic anisotropy systems, there exists a marginal dimensionality to^'^ separating two different 

^Another option would be to consider a random field disorder [ 12 . However, we do not introduce it since an 
appropriate model differs from those considered here by its upper critical dimension and this will make our account 
even broader. The same concerns strong dilution at the percolation threshold. 
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scenarios: for m < mj? dilution causes a change in the universahty class, whereas the asymptotic 
critical exponents of m > mJ?''-component systems remain unchanged under dilution. The so-called 
Harris criterion allows us to determine this marginal dimension from the pure system heat capacity 
behaviour [llTi. It states that critical exponents are not altered by dilution, if the heat capacity of 
the pure system does not diverge (a < 0). From Tabled one concludes that in three dimensions it 
is the Ising model which changes the universality class. 



3.3.2. Random anisotropy 

Another way of introducing randomness to the model H3.1|l is to consider the case, when each 
spin is subjected to a local anisotropy of random orientation. The resulting Hamiltonian reads [ 

n = -lYl J{\^-'R'\)SrSr, - DY^i^RStt)^- (3.6) 

Here, > is the strength of the anisotropy and is a random unit vector pointing in the direc- 
tion of the local anisotropy axis. The random anisotropy model (|3.6|1 is relevant to the description of 
a wide class of disordered magnets. It was first introduced to describe magnetic properties of amor- 
phous alloys of rare-earth compounds with aspherical electron distributions and transition metals 
[ I18|. Today the majority of the amorphous alloys containing rare-earth elements are recognized 
to be random anisotropy magnets [I19|. 

As in the random-site case, the model should be completed by choosing a certain distribution 
p{x-r) for the random variables xr. Most often, two different distributions are considered [ I21|. 
The first is an isotropic one, where the random vector points with equal probability in any 
direction in the m-dimensional hyperspace: 

p(*.,.(/d"«)".^. (3.7) 

Here r(a:) is Euler gamma-function, and the right-hand side presents the volume of the Tri- 
dimensional hypersphere of unit radius. This distribution mimics an amorphous system without any 
preferred direction. The second distribution restricts the vector xr to point with equal probability 
along one of the 2m directions of the axes ki of a hypercubic lattice: 

where S{y) are Dirac (5-functions. This distribution (sometimes called a cubic one) corresponds to 
a situation when an amorphous magnet still remembers the initial cubic lattice structure. 

It is generally believed that ferromagnetism does not exist for the 3d model H3.6|l with the 
isotropic random axis distribution H3.7(l [ 19 . However, anisotropic distribution (|3.8(l leads to a 
magnetically ordered low temperature phase and the transition belongs to the random site Ising 
model universality class [I20|. 



3.4. Frustrations 

An archetype of a model describing the effect of frustrations is the model of stacked triangular 
antiferromagnet [I22|. In this model, the antiferromagnetically interacting spins are placed on the 
sites of 2d triangular lattices stacked in register along the orthogonal direction. The sign of the 
interlayer interaction is unimportant, because there are no frustrations in orthogonal direction. 
The Hamiltonian reads: 

y js^Sn', J^l i'^^' inside a plane, 
2 «• «• ' I J2 , between planes. ^ ^ 

(R,R'> 
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Sum in H3.9|l runs over the nearest neighbours of the above described lattice. Systems which are 
characterized by the Hamiltonian H3.9|l exhibit noncolhnear spin ordering. An example is given by 
the famous 120° structure: each spin in a layer forms 120° angles with the neighbouring spins. 
Although model H3.9|l is formulated for general m, of most interest are values m = 2 and m = 3. 
Namely for these values of m the model has experimental realizations [ .14i i23; and describes 
noncollinear ordering of stacked triangular antiferromagnets as VCI2, VBr2, CsMnBra, and helical 
magnets as Ho, Dy, /?— Mn02. In the latter substances, frustration is induced by the competition of 
ferromagnetic nearest-neighbour and antiferromagnetic next-nearest-neighbour interactions, which 
acts only along one lattice axis: Ji > and J2 < 0, correspondingly. When the ratio J1/J2 exceeds 
a critical value, in the low-temperature phase the spins align ferromagnetically in a plane and 
form a spiral along the orthogonal axis. Moreover, at m = 3 Hamiltonian H3.9(l also describes A/B 
transition in He'^. 

Numerous experimental and MC studies performed so far have not lead to the definitive con- 
clusion about the order of transition into non-coUinear state. There is no unique answer from the 
theoretical viewpoint cither [ 1141 123| . 

4. Renormalization 

In order to explain the main ideas of the renormalization group (RG) theory and to show how 
this method works in practice we shall study the critical behaviour of the Id Ising model by means 
of the RG approach. The explanations given below are due to Michael Fisher who compared 
this study with the Bohr-Sommerfeld picture in quantum mechanics. Passing from the classical 
mechanics to the full account of quantum mechanics, Bohr-Sommerfeld's picture represents only a 
crude approximation. Nevertheless, it introduces some important ideas, just like RG study of the 
Id Ising model enables one to introduce ideas of renormalization and scaling in the critical region. 

4.1. RG transformation 

We start from the Hamiltonian of the Id Ising model in the presence of an external magnetic 
field if, which being divided by fceT is written in the following form: 

N N N 

-n,s^KY,s,s,+i + hY,s, + cJ2'^, (4.1) 

where K = J/ {ksT), h = H/ [k^T)^ and the last term C X^jLi 1 = -^C' is added for the convenience 
of forthcoming calculations. We call this temperature-dependent Hamiltonian an effective one, TLcS- 
Note, however, that contrary to the effective Hamiltonians considered in the forthcoming section|31 
it contains full information about the spin Hamiltonian. Given {K^ h, C} one completely specifies 
Heff' thus it can be regarded as a point in a space of 3 parameters {if, /i, C}. With change of 
{K, h, C}, this point moves. 

One of the first approaches to the RG is to regard it as a specific way of calculating the partition 
function Zjv{H}: 

Si=±l Siv = ±l 

where we have normalized the expression for the partition function by the partition function of 
the ideal model (Z''^<='^^ = 2^) to obtain Zjv{Hcff} ^t^oo 1- The RG method of calculating Zn 
is to step-by-step perform the summation in (|4.2(l in such a way as to try as much as possible to 
preserve the system as it used to be prior the summation. 

Let us do this by means of decimation procedure: performing summation over every second 
spin variable (see figure [21 . To show the changes which are introduced by summation over certain 
spin variable Sq let us represent the total Boltzmann weight in the factorized form: 

e-«c« = . . . ^KS.SoHHS-+So)+c^KSoS^+^h(So+s+)+c p^S-,So)PiSo, S+)--- (4.3) 
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Figure 2. Decimation procedure followed by the spatial rescaling. Summation over every second 
spin variable is performed and then all lengths are rescaled in such a way that the new lattice 
spacing a is equal to the old one a. Now any distance R in the original lattice becomes R' = 1/2R 
in the new one measured in the units of the lattice spacing. 



Summation over Sq which enters only P(S'_,S'o) and P(S'o,S'+) will lead to a new fmiction 
P'{S-,S+), defined by: 

P'(5_,5+) = i P{S-,So)PiSo,S+). (4.4) 

So = ±l 

Factor 1/2 is included in 1)4.4(1 because with each spin eliminated one must remove factor ^ from 
the normalizing factor in H4.2|) . Now the RG idea is to express the new factor P' in the same form 
as the initial one: 

P\S-,S+) = e^'^-^++^'''(^-+^+)+^'. (4.5) 
The new parameters define the renormalized Hamiltonian: 

K^ = K^{K',h',C'}, (4.6) 

and Tigjj has a half of initial spins. The result can be formally written as : 

Kff-7^fa[Hcff], (4.7) 

with a spatial rescaling factor h (in our case 6 = 2, see figure EJ- The resulting number of spins 
N' is connected with the initial one by A''' = N/h, and in the case of d-dimensional system this is 
generalized to: 

N' = N/b'^. (4.8) 

Performing summation (|4.4(l and expressing the result in the renormalized form (|4.5I) one can 
get expressions for the renormalized variables K' ,h', C (it is proposed to the interested reader to 
do this). 

Now the renormalized model is characterized by the Hamiltonian ^{'^^{K' ,h' ,C') having the 
form similar to the initial one, though the lattice constant after taking the sum over every second 
spin is equal to 2a (see figure|21l. In order to have the renormalized model looking as the initial one 
we resettle all lengths in such a way that the new lattice spacing a' is equal to the old one. Now 
any distance R in the original lattice becomes R' = 1/2R in the new one measured in the units of 
the lattice spacing. For the arbitrary b we have the following mapping: 

R^ R' = R/b. (4.9) 
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Now let us have a look at the behaviour of the spin-spin correlation function {SoSr). First let us 
renumber the remaining spins to have the labels arranged in the subsequent order (see figure [2J|: 
S2 ^ S[, 5*4 5*2, S2k ^ S'l^, . . . Ot, taking into account H4.9|l . S2R' = S'^,. Since 

the remaining spins S' have not changed under the renormalization procedure, the renormalized 
correlation function is equal to the original one: 

{S0S2R') = {S'oS'j,,). (4.10) 

It follows that if the original correlation length is <^ = £,{Hcff) then the renormalized correlation 
length is two times smaller. Or, for a general b: 

e(Hoff) = baKs). (4.11) 

The RG procedure has the effect of shrinking the correlation length. Recalling that ^ becomes 
infinite at T = Tc one can state, that the RG procedure is driving a system away from criticality 
(if it was not critical). Already here it is seen that the RG transformation we are considering has 
a deep connection with the critical properties of a system. 



4.2. Fixed points and RG flows 

To proceed in the determination of these properties let us have a closer look at the partition 
function and at the free energy of the system. First, after taking a partial trace of the partition 
function Zi^{Ti.cfi} we get the Boltzmann distribution in N' spin variables S': 

^-KAS'} ^ S'p^';,e-^<=«(^\ (4.12) 

where Spfi„ stands for the trace over TV" = N — N' spin variables S". Now taking trace over the 
remaining N' spin variables one gets: 

Z^4Kff} = ^P^'e-«-(^^ (4.13) 
Substituting into the right hand side of l|4.13|l its explicit form given by (|4.12|1 one gets: 

This relation can be rewritten in terms of the flow equations which describe a motion of a point 
describing an effective Hamiltonian: 

K' = nK{K,h), (4.14) 
h' = nh{K,h), (4.15) 
C' = b'^C + no{K,h), (4.16) 

for the couplings defining HeS- As far as the temperature T enters parameter K, the above relations 
imply the flow equations for T as well. To study this let us consider H — and write the flow 
equation for the temperature as: 

T =>T' = n{T), (4.17) 

with TZ.{T) being appropriate function of T. Suppose that TZ{T) has a form given in the flgureO 
The important feature there is that Tl{T) crosses the line T' = T aX some point T = T* . This point 
is called the fixed point (FP). It is clear that when T is smaller or greater than T* the successive 
application of the renormalization procedure drives the system away from the fixed point. On the 
other hand at T = T* the system remains at the fixed point under the application of the RG 
procedure. Recalling the flow equation for the correlation length H4.11|) which can be written as: 



one has that in the FP this reads: 
Since 6 > 1 it is possible only when 



^(T) = 6e(T'), (4.18) 
^(T*) = b^{T*). (4.19) 
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Figure 3. Temperature renormalization function TZ{T). TZ{T) crosses the line T' = T at the 
fixed point T ^ T* . When T is smaller (T = To) or greater (T = Ti) than T*, the successive 
application of the renormalization procedure drives the system away from the fixed point. On 
the other hand at T = T* the system remains at the fixed point under the application of the 
RG procedure. 



• ^(T*) = oo and 

• ^(T*) = 0. 

The first case characterizes the critical point T* = Tc, whereas the second corresponds to the 
vanishing of the correlation length at zero or infinite temperatures (when the spins are frozen at a 
ground state or totally uncoupled). Note, that this corresponds to two FPs more (T* = 0, T* = oo) 
in the figure |21 

Having identified the fixed point T* of the RG transformation as the critical point of the system 
let us study what knowledge about the critical exponents can be obtained based on the properties 
of the RG transformation. First let us linearize it in the vicinity of the FP introducing the variable: 

. = 1^ = 1^ ,4.20) 

and replacing the plot of TZ{T) near T* by its tangent at T* . 

Then after the renormalization the temperature deviation will be: 

r' = r(i) ~Ai(6)r (4.21) 

for small enough t, where Ai(&) is the slope of the tangent. To find its dependence on b let us 
apply the renormalization procedure twice: 

r" = r(2) ~Ai(&)Ai(&)r. (4.22) 

Such transformation should be equivalent to transforming with a spatial rescaling factor b^: 

t(2) ~ Ai(&2)r, (4.23) 

or 

Ai(6)Ai(6) = Ai(&2). (4.24) 
This leads to the conclusion that Ai(6) should have the following form: 

Ai{b) = b^\ (4.25) 
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with Ai being constant independent of b. Now for the correlation length after / transformations 
one gets: 

e(r)-6'e([Ai(&)]'r) = fe'^(6'^^r). (4.26) 
As far as (|4.2(j|) holds for any / (and for small t(!)) let us choose it to satisfy: 

= 1, 



or 



Then from (|4.26|l one gets: 



6' = (l/r)i/^\ (4.27) 
C(r) = (l/r)i/^^e(l) 



or 



e(T)-r-i/^i. (4.28) 
Comparing (|4.28(l with the definition of the critical exponent i^: ^(t) ~ t^"^ we get: 

j/=l/Ai. (4.29) 

So as we have seen, the knowledge of the linearized RG transformation enables one to determine 
the critical exponent! 

Note, that having carried out similar analysis for the free energy per spin. / = — A^^^ln^Ar, 
one arrives at the expression: /(r) = 6~'^'/(6^i't). Again, choosing b from H4.27(l one gets /(r) = 
T'^/^^f(\) = t'^^ f{l). Comparing the latter with the scaling behaviour of the free energy in the 
vicinity of the critical point /(t) ~ r^^" one proves the hyper scaling relation: 

2-a^dv. (4.30) 

Similarly, considering the case of non-zero magnetic field one can write the recursion relations 
for the renormalized temperature and field. Their linearization in the vicinity of the fixed point 
leads to the matrix of the linear RG operator with two eigenvalues, defining two different critical 
exponents. 

To complete the RG study of the Id Ising model it remains to perform the above described 
procedure explicitly and to find the values of the critical exponents (at zero-temperature FP T = 
H — 0, [1^). The above described transformations concisely give the main features of the RG 
transformation, allowing us, in particular, to define the critical exponents of the system. 



5. From the spin Hamiltonians to the effective ones 

Very often a starting point for the RG study of critical behaviour of a many-particle system 
is its effective Hamiltonian. Taken, that in the model description the system is identified with the 
spin Hamiltonian (see section , the effective Hamiltonian arises as a certain "metamodel" . It 
shares global features of different spin Hamiltonians: their dimensionalities, symmetries, type of 
interparticle interaction and as a result, brings about common features of their critical behaviour. 
Below, we shall show how to obtain the effective Hamiltonians for the spin models of section |31 



5.1. m-vector model 

To get an effective Hamiltonian we proceed as follows. Let us define the free energy !F and the 
partition function Z of a spin model H3.1|l as: 

T=-fi-HnZ, Z = Spc-'3", (5.1) 

where /? = (fceT)^^ and Sp(. . . ) as usually means the sum over all possible states. In our case it 
corresponds to the integration over the surface of m-dimensional hypersphere (we take it to be of 
unit radius): 

Sp(...)=n / dSRdi\Sn\-l)i...). (5.2) 

R 
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Let us introduce the Fourier-transforms of the variables 5r by: 

S^-^E^'"'"^^' ^^-4^Y.^''''''Sn. (5.3) 



Here and below, when it will not be mentioned explicitly, vector R spans all sites of the lattice 
whereas k changes in the first Brillouin zone. Now the Hamiltonian (|3.1f) can be rewritten as : 

W = -i^i^(fc)^k^_k, (5.4) 

k 

where we have introduced Fourier transform of the potential of interaction by: 

J^^^ = ^ E '^(^) = E c-'''''^(i?)- (5.5) 

k R 



Now the partition function (|5.1|) reads: 



/3 



In order to take the trace in H5.6|) let us transform it to the expression with linear dependence on 
Su_. To this end, one makes use of the Stratonovich-Hubbard transformation introducing the field 
variable 0k which is conjugated to the spin variable Sk by the identity: 

e-P^ = JJe^^^--^^-- =]J( ^^^^^^^^ )'"/^ J d4e^^^-^^^-+^^''^^-- 



k k 

'"(d0)e^''^WTl'?''l'+^'"?->'). (5.7) 

As far as the Stratonovich-Hubbard transformation was performed for each function 5k we arrived 
at the functional integral J{d$) = nRn"=i I^oo "^"^R '-'^^^ &e\d variables 



R 

From now on we omit the coefficients in front of the functional integral. 

We have reached our goal: now the trace in the expression for the partition function concerns 
only the last term of the exponents under integration in (|5.6(l which is a linear function of S: 

Z = Spe^^k-^e^')^"^-" - y"(d0)e^^K /J^^I'^-l'Spe^R^^^'^^. (5.9) 

The last step is to take integral over 5r, (recall that the trace is defined by equation (|5.2|l ). This 
can be achieved by passing from the m-dimensional Cartesian coordinates , , . . . , to the 
m-dimensional polar ones S,6i, . . . , Om-i- Performing this integration and representing the result 
as a series in one gets for the partition function: 

(d0)e^^''(:^"""^l'^''l'"^i^^'=2 (Wl-^J^I", (5.10) 

where the coefficients U21 readily follow: U2 = — 1/m, W4 = 6/nr^{'m + 2), ... Expression (|5.10(l 
gives the functional representation of the partition function of the TO-vcctor model. So far we have 
not gained a lot by the above described transformations: the difficulty of taking Sp of the initial 
expression (|5.4|l now is transformed into the difficulty of calculating a functional integral H5.10|l . 
However, in order to study critical behaviour of the model (|3.1|) . expression (|5.10l) can be further 
approximated. It appears that all of higher than the fourth powers do not effect the asymptotic 
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critical behaviour at d = 3: they do not change the value of the fixed point of RG transformation 
(cf. section 0} and are irrelevant in the RG sense [3,4. Therefore one can be restricted to the 
0^ model. Being interested in the long-range correlations arising in the system in the vicinity of a 
critical point, one substitutes the Fourier image of the interaction potential by its expansion for 
small wave vector values k: v{k) ~ 1^(0) — 1 / 2\v" {Q)\k'^ and writes the contribution fc^(/)k'/>-k in 
the form X]r(^'/'R.)^ re-scaling variables to get the gradient term without any coefficient. One 
further passes to the continuous limit — s- J dR and gets for the free energy: 

~ In j (d0)e-"-ff (5.11) 

with the effective Hamiltonian: 

Heff = |d^i?|i((V0)2+/.g0^) + ^</-4|, (5.12) 

where </> = 0r, 0^ = |(/)rP, variables fiQ and uq are called bare mass and coupling. From the above 
derivation we know that uq is positive (uq ^ from equation (|5.1U|) ). which ensures the existence 
of the integral H5.11|l . Effective Hamiltonian (|5.12() shares common global properties with the spin 
Hamiltonian H3.1|l : the dimension of space (i.e. dimension of vectors k, R), dimension of the order 
parameter (m) and its symmetry: the functional representation we obtained is symmetric under 
group of rotations in the m-dimensional space 0(m): it depends only on \(f)\. Let us now see what 
differences will appear in the functional representation of the non-ideal models. 



5.2. Cubic model 

It is intuitively clear that an effective Hamiltonian of the cubic model H3.2|) will differ from that 
of the m-vector model (|5.12|) by the presence of terms with new, cubic symmetry. Indeed this is the 
case that one can easily check following the next steps: starting from the Hamiltonian (|3.2|) via the 
Stratonovich-Hubbard transformation one obtains the functional representation for the interaction 
part (which coincides with equation H5.7|l 'l. An additional term e~^^^'^ can be represented in 
a form of the functional derivative e~^^^i'^ Z^*-"^ -* , resulting in the following expression for the 
partition function of the cubic model (cf. equation l|5.9|) '): 

Z - y (d<^)e'^^'':5T;T^)l'?''l%-/5^i:i^i:.9V9(0k)*Spe^i^^^'^=^. (5.13) 

Taking trace in equation (|5.13|) leads to the familiar expression H5.10|l . However, the derivative 
9^/9(01^)* gives rise to the contributions of cubic symmetry, the lowest order contribution being 
proportional to voJ2i{4'\i)'^ ■ The resulting effective Hamiltonian reads: 

n.^ = |d'^i?|l((V0)2+^^^2)+^<^4 + ^g(</>^)4|, (5.14) 

and contains two bare couplings uq and vq, corresponding to the (f>'^ terms of different symme- 
tries. Coupling Uq is positive, whereas the sign of the coupling vq coincides with that of V in 
equation (|3.2|) . At = one recovers an effective Hamiltonian of the m-vector model (|5.12ll . 



5.3. Weakly diluted quenched m-vector model 

The peculiarities of the free energy calculation for the model (|3.3|l consist in averaging over 
quenched disorder. Indeed, for each configuration of empty and occupied sites in (|3.3|l one can 
write a corresponding configuration-dependent partition function Zconf'- 

2conf - Sp^e-^^, (5.15) 
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where Spg concerns spin degrees of freedom and is defined by H5.2|l . The free energy is obtained as 
the configurational average: 

^=-/5~l(lnZeo„f)conf. (5.16) 

One of the ways to proceed is to make use of the repUca trick [ [53] , which aUows us to avoid 
integration of the logarithm in 1)5. 16() substituting it by a power function: 

Z" — 1 

lnZ=Um . (5.17) 

n^O n 

Then Z'^^^^ can be written as: 

Z,"„„f = Spe^^i^.t^"^(l^-^'l)^"i^^'^^'. (5.18) 

with obvious notations Sp(. . . ) — Ha^i Sp^Ja (. . . ) and = crS'^^. Upcoming calculations closely 
follow lines of the subsection 15.11 Introducing by the Stratonovich-Hubbard transformation field 
variables cj)^, conjugated to one is able to take the trace over spin subsystem and is left with 
the configuration-dependent partition function: 

Note that cr appears in (|5.19)) in the first power because from cr = {0, 1} it follows that (cr)' = 
cr. The last step is to perform configurational averaging of 1)5.19)1 with the distribution function 
()3.5)l . Let us represent the result of averaging in the exponential form 



conf -'■ -'- 
R 

where Kp(c) are cumulants of random variables cr and can be easily calculated for the random 
variable distribution l),'15)l : 

CR = {0,1} 

«2 = (cr)Vcr)-( crp(cr) =c(1-c),.... (5.21) 

CR = {0,1} \CR = {0,1} / 

For the free energy one gets: 

T = -/3-i(lnZeo„f)conf = i™o {tJ" y"(d0)"e-^['^l - l/n| , (5.22) 

where the (n-dependent) coefficient ci can be recast explicitly following all the steps of calculations 
described above and the free energy functional F[(j)\ is given by: 

k a=l ^ ^ ' pjsl ^ R = l ^ ' a = \ / 

Again, as in the previous subsection, restricting ourselves to the 0'* terms, expanding the short- 
range interaction potential, re-scaling the fields and passing to the continuous limit one gets for 
the free energy: 

T r-. \ (d0)e~^<=« , (5.24) 



14 



Introduction to renormalization 



now the proportionality sign also hides, besides the coefficient, the replica limit, as written explicitly 
in (15.2211. The effective Hamiltonian reads: 



/[ ri n n | 

d'^^UE[(^'^")'+/^o(m + 9E(^")' + 5 E im<P'r\- (5.25) 

[ Q=l ■ a=l ■ a,/3=l J 

Here, the coupling uq is positive (being proportional to cu4, 15.f 0|l ') whereas the coupling vq is 
proportional to c(c — l)u2 and is negative. The last term in (|5.25l) is present only for non-zero 
dilution: it is directly responsible for the effective interaction between replicas due to the presence 
of impurities. 

5.4. Random anisotropy model 

To treat the random anisotropy in the spin Hamiltonian , one first writes the configuration- 
dependent partition function for a fixed local anisotropy axes configuration {x}: 

2conf({i}) - Sp^e-'^". (5.26) 

As in equation H5.15(l , trace in equation H5.26|l concerns only the spin degrees of freedom. Applying 
the Stratonovich-Hubbard transformation to the interaction part of the Hamiltonian one presents 
(|5.26|l in the form of the functional integral and is able to take the trace: 

2conf({x})^ j{d$)e-^^^\ (5.27) 

with 

n{x) = j d'^i?|i [{V(j)f + filq^^] - Di{<j)x f + vicb^ + 2i02(0i)2 + . . . |. (5.28) 

Again, the expansion for small k was performed and the continuous limit has been taken. Explicit 
values for the coefficients fii, Di, vi, z\ are given in [ I2(J) . 

The rest of calculations follow the steps outlined in the section 15.31 for the quenched disorder 
the free energy is defined by (|5.1()l) , where averaging is performed over the random axis distribution 
(given by (|3.7() or H3.8|l ). Subsequently, one substitutes In Z by a power function via the replica 
trick (|5.17(l . For the isotropic random axis distribution H3.7|l one gets the effective Hamiltonian [ 
12 1| containing three couplings of different symmetry: 



/[ n n m 

[ a=l a,p=\i,j=\ 



Here and below ip^ = J2a=ii't'°')'^ ■ One can check the signs of the couphngs: uq ~ m^(m+2) ^ ^' 
Vq ^ D > 0, wq ^ ^(^^_2) < 0- Moreover, from the explicit expressions for the couplings Wq and Uq 
one gets for their ratio Wq/uq ~ — m. The latter relation determines a region of physically allowed 
initial values in the (u — v ~ w)-space of couplings. 

For the cubic distribution H3.8|l the average over the random variables {x} leads to the effective 
Hamiltonian [ I21j with four couplings: 



'HcS — 



a=l 

2 ™ " 



m n 



+W0 J2 E (-^n' (-^f ) + yo E E (-^n' } ■ (5.30) 

i—l Q,/3— 1 i—l ct—1 
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The signs of the coupHngs are as follows: uq ^ D"^ /{2rn?) > 0, ~ -D > 0, ~ — D^/(2m) < 0. 
The last term in H5.3U|I is of cubic symmetry. It has to be included since it is generated if one further 
applies the RG transformation. Therefore j/o can be of either sign. The symmetries of wq terms in 
1)5. 29(1 and H5.30() differ. However the ratio wq/uq = —to holds for the effective Hamiltonian (|5.30() 
as well. 

5.5. Stacked triangular antiferromagnet 

It is an interaction potential H3.9|l that makes a difference in a derivation of the effective 
Hamiltonian for the frustrated magnet and the regular one. One can find the details of the procedure 
to derive an effective Hamiltonian via the Stratonovich-Hubbard transformation in [E2. The main 
difference from the above described calculations for the to- vector magnet f subsection 15.111 arises 
when one proceeds with the Gaussian term in equation (|5.1()(l : 

5 "0 " ? ^'-^'^ 

The sum over k in l|5.31|l spans the 1st Brillouin zone. For the to- vector model, d(k) has a minimum 
at A: = and an expansion in k around this minimum results in the effective Hamiltonian H5.12|l . 
Note, that one minimum in d{k) leads to one field variable 0k in the effective Hamiltonian. For 
the frustrated model, the d{k) has a maximum at fc = and two independent minima at ^ in 
the 1st Brillouin zone. Subsequently, the zone can be rearranged into separate subzones and the 
field 4>k over the original Brillouin zone can be decomposed into two fields, (/)i_k and 02, k, with k 
confined to the new subzone. Now, the minima of d{k) occur at the origin and pertain to the fields 
01 and 02 ■ The above decomposition results in the appearance of the terms of new symmetry in 
the effective Hamiltonian, which now reads: 

Heff = / d"^^ [m^(0? + 4) + (V0l)^ + (V02)^] + ^ [0? 0^]% ^ [(0, . 0,)2 - fyj] I , 

(5.3 2) 

with a scalar product of TO-component fields 0i ■ 02 = '^"Li The coupling uq in (|5.32() is 

positive, whereas the sign of the coupling vq determines a type of spin ordering: a non-coUinear 
order occurs for vq > 0. For vq < the fields 0i and 02 tend to allign parallel or antiparallel, this 
corresponds to the sinusoidal ordering or the linearly-polarized spin-density wave. 

6. RG explanation of criticality in non-ideal systems 

Now with the effective Hamiltonians (section |5)l and the RG method (section at hand one 
can proceed further, applying the method in order to study the critical behaviour of the models. 
However, to be honest with the reader, he should be warned that the whole story is not that simple 
as shown in the section 01 It is not only because the very models we are interested in are much 
more complicated as compared to the Id Ising model: in what follows below we will be interested 
in the criticality of 3d systems. The real space renormalization described in the section 0] works for 
the low-dimensional systems, whereas for realistic 3d systems one rather applies the RG equations 
in the momentum space. Moreover, the very approach to RG as to the way of calculating the 
partition function of the system [ 0] has its alternative known as the field-theoretical RG [ 0] . 
The latter is a tool to cope with the divergencies of correlation functions (vertex functions). Here, 
the RG procedure consists in the controlled rearrangement of the perturbation theory expansions 
giving rise to the RG equations. However, the underlying notions of both procedures are the same: 
given the effective Hamiltonian one applies the RG transformation and studies the fiow equations 
for the couplings of the Hamiltonian. If the fiow equations possess the fixed point (FP), which is 
stable and reachable from the initial conditions, it corresponds to the critical point of the system. 
Scaling arises in the vicinity of this FP and the universal values of the exponents governing scaling 
of different physical quantities may be found. Having this preamble in mind, we can make a brief 
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review of what resuhs for the RG flows, FPs and exponents of the non-ideal models have been 
found so far and how they differ from those of the "ideal" to- vector model. 

6.1 . m-vector magnet 

Put in a formal way, the question: does the 3d m-vector magnet exhibit a critical point (a Curie 
point for a ferromagnet or a Neel point for an antiferromagnet), transforms into the question: is 
there a stable and reachable FP in the flow equations (cf. H4.14ll - H4.16|l 'l for the couplings of its 
effective Hamiltonian? The flow equation for the coupling u of the effective Hamiltonian (|5.12|l can 
be written in a general form of an ordinary flrst order differential equation: 

n 7/ 

f3{u), (6.1) 



dln£ 



with the RG flow parameter £ and P -function P(u). The parameter £ may serve to measure a 
distance to the critical point: ^ corresponds to T — > Tc, and a specific form of the /3-function 
depends on a choice of the RG procedure. The data discussed in this chapter have been obtained 
within the field-theoretic RG [|4j approach at d = 3. Giving numbers, we shall refer to the source, 
where they were obtained. However we shall try to be not too specific, mentioning several review 
papers for the interested reader. 




(a) (b) 
Figure 4. /3-function of the 3d m = 1 model in successive perturbation theory orders ranging 
from 1 to 5 as shown by the labels in the figures. Left: naive evaluation of the function 16.211 . 
Right: resummation of (lfci.2ll taking into account asymptotic properties of the series. Note a 
decrease of the difference between FP coordinates found in successive orders of the resummed 
perturbation expansion, figure [IJs. 

To give an idea about the RG expansions and their treatment, we write down several first terms 
of the /^-function (EIJ [EB]: 

P{u) = -u{e -u + 3u^(3to + 14)/(to + 8) + • • • ), (6.2) 

with e = 4 — d. For a quantitative analysis, one can develop an e- expansion looking for solutions 
of the FP equation 

l3{u*) (6.3) 

in a form of a series u* — ^^u'^'^e'. Alternatively, one can solve the non-linear equation H6.3I) 
directly at d = 3 (e = 1). However, as it is well known by now, the perturbative RG expansions 
have zero radius of convergence and are asymptotic at best (cf. behaviour of the function (|6.2|l . 
figure 2^). Special procedures of resummation have been elaborated to deal with them. We give 
an example of how one of these procedures works transforming a divergent series, figure 0^, into 
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a convergent one, figure^jD. In the last figure, the function H6.2(l has been resummed by the Pade- 
Borel resummation. The procedure consists of several steps. First, assuming the factorial growth 
of the coefficients Ci in the expansion P{u) = CjW* H6.2() . one constructs the Borcl trransform of 
the initial function /3(u) via: 

/^^H = E^- (6-4) 

i 

Then, the Borcl transform is extrapolated by a Pade approximant [K/L]{u). The last is the ratio 
of two polynomials of order K and L such that its truncated Taylor expansion is equal to /3^(u). 
The resummed function is then calculated by an inverse Borel transform of this approximant: 

/•oo 

P{u) = / dtexp{-t)[K/L]{ut). (6.5) 
Jo 

Similar techniques are currently widely used in analyzing the RG expansions [I^I27|. In particular, 
the numerical estimates of different physical quantities given in this chapter (and in Table ^ as 
well) have been obtained using resummation techniques. 

One sees the presence of two FPs in figure 0Jd: an unstable Gaussian FP G u* = (the 
slope of the /^-function is negative, df3{u)/du\u* < 0) and a stable one, u* ^ 0. Here, the slope 
of the /9-function is positive: starting form any initial conditions with m > the solution of the 
differential equation (|6.1(l reaches its FP value. This FP corresponds to the critical point Tc of 
the 3d Ising model. Similar behaviour of the /3-function is found for other values of m, therefore 
the FP with m ^ for general m is called a 0(m)-symmetric or Heisenberg FP H. We do not 
show the procedure of calculating the critical exponents: the other RG functions being evaluated 
at this FP bring about the asymptotic values of the critical exponents (in particular, those given 
in Table ^ . These exponents govern criticality of systems of different microscopic nature, which 
can be described by the effective Hamiltonian (|5.12|l . It is said that these systems belong to the 
0{m) universality class. In the RG picture, different microscopic origin is reflected in different 
initial conditions for the flow. However, the FP location and stability is deflned solely by the global 
features: dimensionality, symmetry, interaction type.'^ Systems which share the global features 
belong to the same universality class. 



6.2. Cubic model 

The main difference of the effective Hamiltonian of the cubic model (|5.14() as compared to 
that of the m-vector model (|5.12|) is that it contains one more coupling v of different symmetry. 
Therefore, two /3-functions describe the RG flow: 

Puiu,v), ^ f3^(u,v). (6.6) 



dln£ ' ^ ' dln^ 

Stability of a FP is now defined by the stability matrix composed of the /3-function derivatives 
dPu,v/d{u,v}. The FP picture, which arises form an analysis of the /3-functions at d = 3 [|^| is 
schematically shown in figure [3 Four FPs are obtained: unstable Gaussian G {u = v ^ 0) and 
Ising I (u = 0, u 7^ 0) as well as Heisenberg H {u ^ 0,v = 0) and mixed M {u ^ Q,v ^ 0). The 
stability of the two last FPs depends on the value of m. 

Two different regimes for the RG flows are observed. For small m < m^^^ the FP H is stable. At 
this FP the system does not feel a presence of the cubic coupling, v = 0. Therefore, for m < m^"*^ 
the cubic model belongs to the 0{m) universality class: its exponents coincide with the exponents 
of the m-vector model, table ^ However, with an increase of m the FP M approaches H and 
at m = m^^^ both FPs coincide: a crossover to the new regime occurs. For m > m™'^ the FP 
M becomes stable and governs the critical properties of the cubic model in the new universality 
class. The marginal dimension value is slightly less than three: m™'' — 2.862(5) [|2H|- From this 



^In particular, here we consider systems with the short-range interaction. 
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Figure 5. FPs and RG flows of the cubic model. Unstable FPs are shown by discs, stable FPs 
are shown by squares, (a): m < rri^^ , the Heisenberg FP H is stable, (b): m > nri^i^^ , the mixed 
FP M is stable. The dotted line in figure |^ shows a separatrix: the RG flows which start below 
this line do not reach H. Note that for m > nf^^ the u-axis is a separatrix: all flows that start 
from the initial conditions with u < do not reach any FP (run-away solutions). 

estimate it follows in particular, that a cubic Heisenberg (m = 3) magnet does not belong to the 
0(3) universality class. Its critical exponent being estimated in the fixed point M read [I29j: 



7 = 1.390(12), J/ = 0.706(6), r; = 0.0333(26), /3 = 0.364(15), a = -0.118(18). (6.7) 



Numerically, these values are close to their counterparts for the m- vector model (cf. Table ^ . 
However the principal difference arises from the above analysis: as one can easily check solving 
the system of differential equations (|6.6|l for m > m™*^ the RG flows with w < cannot reach the 
stable FP (an abscissa t; = serves as a separatrix for the flows). As it follows from sections l3.2ll5.2l 
negative v corresponds to ordering along diagonals of m-dimensional hypercube. Therefore, the RG 
analysis results in a statement that such ordering cannot occur via a 2nd order phase transition: 
ferromagnetic crystals with three easy axes should undergo a 1st order phase transition. It is 
worth noting here, that whereas the presence of stable and reachable FP brings about the 2nd 
order phase transition, its absence signalling only that the 2nd order phase transition does not 
occur. The nature of the low-temperature phase and the scenario of how it is attained remains to 
be checked by other methods. 

6.3. Random-site dilution 

Formally, an analysis of the weakly diluted quenched m- vector model resembles those we dis- 
cussed in the former subsection 16. 21 Indeed, both effective Hamiltonians H5.14|l . (|5.25|) contain two 
couplings of different symmetry and give rise to the already familiar FP picture. Moreover, the 
effective Hamiltonian (|5.25f) at m = 1 coincides with H5.14|l .'^ However, our goal is to analyze it in 
the replica limit n = 0. Note, that now the physically meaningful values of couplings are m > 0, 
v < 0, see section l57^ The FP picture and the RG flows at d = 3 are shown in the figure El [ il3j . 
Besides the familiar FPs G, H, and M a polymer FP P is present. It is stable and corresponds 
to the 0(m = 0) universality class, however it is never reached from the initial conditions w > 0, 
u < 0. 

Again, the new marginal dimension mj?'' governs the crossover between the new and the 0{m) 
universality classes: FP M is stable for m < mj?''. A search for the value of mj?'' shows, that at 
m = mj?'' the heat capacity critical exponent of an undiluted system changes its sign: a{m^^^) — 0. 
In this way one recovers the Harris criterion section l3.3.1l translated into the RG "language". 

*One can check it by further substitution n — » m, {uo,vo} {vo,uo}. 
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(a) (b) 
Figure 6. FPs and RG flows of the weakly diluted quenched m-vector model. Unstable FPs are 
shown by discs, stable FPs are shown by squares, (a): m < mj?'', the mixed FP M is stable, (b): 
m > rric'', the Heisenberg FP H is stable. Flows 1, 2, 3 of the figure|SJj are further treated in 
figure Eb- 

The numerical value of m^'^ being slightly less than two, m^'' — 1.912(4) ["30", only the Ising model 
(m = 1) changes its exponents upon dilution. Indeed, the numerical values of the exponents read [ 
EH: 

7 = 1.330(17), 1^ = 0.678(10), = 0.030(3), /? = 0.349(5), a = -0.034(30) (6.8) 

and differ essentially from those of the 3d Ising model (table 

However, both in the experiments and in the MC simulations one deals with the system not yet 
in an asymptotic region, where exponents do not attain their FP values and the effective exponents 
are found. Being non-universal, they can be calculated in the RG treatment as functions of the 
flow-dependent couplings. In figure[3 we show an effective critical exponent measured recently 
for the ac susceptibility of the a-Fe86Mn4Zrio amorphous alloy and compare them with the 
theoretical RG calculations of the effective exponents of the weakly diluted quenched m = 3 
model. The exponent was calculated along different RG flows labeled by numbers in figure 
This calculation serves an example how the non-asymptotic effects may be taken into account in 
the RG analysis. Although direct correspondence between the temperature distance to the critical 
point T and the RG flow parameter i is problematic, the RG serves as a useful tool of accompanying 
the studies of effective critical behaviour. 

6.4. Random anisotropy 

On this example we shall show how the RG predicts two different phenomena occuring in the 
random anisotropy magnets. It appears that the type of local random axis distribution crucially 
effects an origin of the low-temperature phase in random anisotropy systems [ I19| . First, we consider 
the results obtained for an isotropic distribution H3.7|l . which leads to the effective Hamiltonian 
H5.29|l [ 1^ . It contains three coupHngs, u, v, w, hence three /3-functions define the RG flows. Solving 
the FP equations one arrives at the FP picture shown in figure |Ht- [ 1211 E]] . 

Let us recall (cf. section I5.4|l . that physically meaningful values of couplings are u > 0, 
V > 0,w < 0. Therefore only the FPs located in this region are shown in the figure. However, 
there is another condition found for the ratio of couplings: wq/uq = —m. The region of typical 
initial conditions to study the RG flow is shown in figure |H1 by a cross. The only stable FP found, 
a polymer FP P is not reachable from the initial conditions. The run-away solutions of the RG 
equations bring about an absence of a 2nd order phase transition. 

A different picture is obtained for a cubic random axis distribution (|3.8|l [ I20| . Here, the effective 
Hamiltonian (|5.30l) contains four couplings of different symmetry, u, v, w, y. The physical initial 
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Figure 7. Effective critical exponent for the magnetic susceptibility of the weakly diluted 
quenched Heisenberg magnet, (a): as a function of the distance to the Curie point r = (T—Tc) /Tc 
in experimental measurements for a-Feg6Mn4Zrio amorphous alloy [ I82| . (b): as a function of the 
RG flow parameter in theoretical RG calculations for the diluted m — 3 model [1151. Different 
curves correspond to different amount of disorder. Note that in the asymptotics (r — * or £ ^ 0) 
the exponent attains its universal value. 

values for the couplings lay in the region fsection 15.4(1 : u > 0, u > 0, w < 0, w/u = —m. A typical 
FP picture is shown in figure IHb for v — 0. Similar to the former case of isotropic random axis 
distribution, the stable FP P cannot be reached. However, one more stable FP M is present. It is 
reachable for the RG flow that starts from the initial conditions marked by a cross in the figure. 
This FP is a FP of the random-site Ising model (subsection 16.311 for any value of m. It means that 
the ferromagnetic 2nd order phase transition in the m- vector magnet with the cubic random axis 
distribution belongs to the universality class of the random-site Ising model and is governed in 
asymptotics by the exponents 1(6. 8|l . 

6.5. Stacked triangular antiferromagnet 

Again, as in the former subsections 16. IhlH^ the RG answer about a possibility of a 2nd order 
phase transition in the 3d stacked triangular antiferromagnet would be a presence of a stable 
accessible FP for the couplings u,v of the effective Hamiltonian 1(5.32(1 . The model possesses a 
rather complicated FP structure, sketched in figure El 

FP picture changes with m and one finds three marginal dimensions mf™ that govern its 
topology. For large m > mf^" the stable chiral FP C+ is present and it can be reached from the 
initial values of the couplings u,v ^ 0. Therefore, the phase transition to the non-collinear chiral 
state is of the second order. At m = mf^", C+ merges with the anti-chiral FP C_ and disappears, 
only one unstable Heisenberg FP H is found for mf^" > m > mf^". With further decrease of m, 
two more FPs appear at m = mf^", S_ and S+. The last one is stable, however it describes the 
sinusoidal phase, which occurs for v < 0. Finally, at m = mf^" this FP merges with H and looses 
its stability with further increase of m. It follows from this analysis, that the 2nd order phase 
transition into non-collinear phase can occur only in the magnets with m > mf^" . 

In principle, the above described FP picture has been known since the first RG studies of 
the problem [ I22( . However, reliable numerical estimates for marginal dimensions mf^" have been 
obtained only recently. The question of interest is, what is the order of the phase transition at m = 
2; 3 when the model has its physical realization? The estimates: mf^" — 6.23(21), mf^" — 1.99(4), 
^chir _ x_43(2) [[23 clearly rule out the possibility of a 2nd order phase transition for m = 2; 3: 
the FP picture is shown in the second figure of figure [HI and no accessible FP is found. 

Note, however, certain controversy in the RG studies performed so far: whereas the above FP 
picture is confirmed by the perturbative RG expansions and the non-perturbative RG approach [ 
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Figure 8. Fixed points of the random anisotropy model. The filled boxes show the stable fixed 
points, the cross denotes typical initial values of couplings, (a): isotropic local random axis 
distribution. The stable FP P cannot be reached by the RG fiow which starts from the region 
shown by a cross in the figure, (b): cubic local random axis distribution. The stable FP P cannot 
be reached, but the random Ising FP M is stable and reachable for the RG flow. 

1231 , an analysis of the FP equations by direct solution of the non-linear FP equations for the 
resummed /3- functions brings about a presence of the stable reachable FP for m = 2; 3 The 
last is associated with the critical point of the 2nd order phase transition. 

7. Conclusions and outlook 

How do the changes in structure effect the critical behaviour of the matter? We tried to give 
an answer to this question taking a 3d m- vector model as an ideal reference system and showing 
what will happen to it under an effect of different non-idealities. The examples considered include 
anisotropy, structural disorder, frustrations: the features one often encounters dealing with realis- 
tic condensed matter objects. The response of a system to such non-idealities appears to be very 
different, ranging from an insensitivity (the cubic m — 2 magnet remains in the 0(2) universality 
class), through softening (the heat capacity of the random site m = 1 magnet does not diverge) 
to disappearance of a 2nd order phase transition (isotropically distributed local random axis de- 
stroys long-range order). These various asymptotic features are accompanied by a complicated 
non-asymptotic effective critical behaviour. 

It is astonishing that all this bunch of phenomena can be explained and accurately described 
within one theoretical framework, the RG approach. An application of basic RG notions of flows, 
fixed points and their accessibility, marginal dimensions and crossovers supported by an elaborate 
machinery to perform and analyze the RG transformation resulted in a coherent picture of phe- 
nomena in the vicinity of a critical point. A lot remains to be cleared up in this picture. Maybe a 
participant of the school or a reader of these lectures will decide to make his or her contribution? 
Good luck! 
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Figure 9. FPs and RG flows of the stacked triangular antiferromagnet model. Unstable FPs 
are shown by discs, stable FPs are shown by squares. Three marginal dimensions mi*^" , m!^" ^ 
7713'""^ govern the FP picture. 
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Questions and answers 

Q (Alexander Lopez): How do the signs of the couplings u and v in the effective Hamiltonian 
H5.14|) determine a type of the low-temperature ordering? 

A Neglecting fluctuations (taking a function (j){r) to be just a variable (p) you can think about 
the effective Hamiltonian as of the Landau free energy. Now, let the reference system display 
a 2nd order phase transition. This means that u > 0. Minimizing Landau free energy and 
looking for the spontaneous magnetization one finds at T < Tc two different non-trivial 
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solutions: = {(ji/^/m, . . . , (j)/y/rn), it exists at v < 0, and = (0, 0, . . . , 0) at u > 0. They 
correspond to two types of ordering: along the diagonals or along the edges of a m-dimensional 
hypercube. For to = 3, these are directions [111] and [100] correspondingly. 

Q (Bertrand Berche): What are experimental realizations of the cubic model? 

A I have already mentioned ferromagnetic crystals. Besides, at m = 3 the model describes a 
ferroelectric phase transition which occurs in SrTiOa at 105 K ( Cowley R.A., Bruce A.D., 
J. Phys. C, 1973, 6, L191). Moreover, since the model provides an example of a system with 
an arbitrary weak first-order phase transition, it is also used as a testing ground to describe 
an elecroweak transition in the early Universe (P. Arnold, S Sharpe, L. Yaffe, Y. Zhang, Phys. 
Rev. Lett., 1997, 78, 2062). 

Q (Carlos Vdsquez): In 1983, Weinrib and Halperin proposed a model to describe an effect of 
extended (correlated) disorder on magnetic 2nd order phase transition. There, the impurity- 
impurity correlation function decays for large separations as g{r) ~ r"". Currently, there exist 
two different predictions for the critical exponents of such a model. What are the methods 
used to obtain them? 

A Indeed, the original result of Weinrib and Halperin ( A. Weinrib, Halperin B.I., Phys. Rev. B, 
1983, 27, 417) was obtained in the first order of the expansion in e = 4 — d, 5 = 4 — a. Recent 
two-loop estimates (V. V. Prudnikov, Prudnikov P.V., Fedorenko A. A., J. Phys. A, 1999, 32, 
L399) are due to the fixed d, a RG technique. Qualitative answer of both approaches is that 
for a < d the disorder is relevant if the correlation length critical exponent of the system 
without defects obeys v < 2/ a. 

Q (Bertrand Berche): Could you comment on logarithmic corrections to the scaling laws? Is the 
critical behaviour at marginal dimensions you were speaking about (e.g. at mj?'' ~ 1.91 for 
the 3d diluted TO-vector magnet) governed by such corrections? 

A Logarithmic corrections arise at upper critical dimension d^^: that is, at the space dimension, 
above which the mean-field theory holds^. For the models I was speaking about, the upper 
critical dimension is four. Indeed, the divergence of an isothermal susceptibility of the TO- 
vector model at d = 4 is governed by the mean field exponent 7 = 1. However, the power 
law singularity is accompanied by a logarithmic one: 

X-i^|T|-iln|T|-("+')/("+«\ (7.1) 

as first derived in: Larkin A.I., Khmelnitskii D.E., JETP, 1969, 29, 1123. In the RG scheme, 
one can see the origin of such corrections solving the flow equation (|6.1|l at d = 4 {e — 0): 

(7.2) 



dln^ 



here the right-hand side is the leading term of the /3-function (|6.2f) . The solution of equa- 
tion ifT?^ 

u = + const, e^O (7.3) 

\\n£\ 

being substituted into an appropriate expression for the susceptibility leads to the above 
dependence H7.1|l . Now, let us return to the marginal dimensions TOc we were discussing in 
these lectures. They are field dimensions and we estimated them at space dimension d = 3. 
As far as the space dimension was lower than d^^ = 4, the logarithmic corrections do not 
appear. Returning back to the example given by formulas (|7.2|) . 17.3|l : even at m — TOc, the 
first power of couplings is present in the right-hand side of functions l|7.2|) for d = 3 and the 
solutions will rather behave as m ~ £const^ 



^Note added in proof: We do not discuss here the logarithmic corrections appearing at low dimensions as those 
in 2d diluted Ising model or 2d g = 4 Potts model. For a list of systems where logarithmic corrections appear and 
for the scaling relations between them see: Kenna R., Johnston D.A., Janke W., Phys. Rev. Lett., 2006, 96, 115701. 
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Q (Dragi Karevski): Do you know other citeria, similar to the Harris one, which predict changes 
in the critical behaviour caused by different types of disorder? For example, what happens 
when disorder is coupled to the order parameter? 

A Indeed, Harris criterion concerns the systems, where disorder is coupled to the energy density 
(look for example at the Hamiltonian (|3.3|l . where random variables cr are coupled to the 
product of spins). It states that critical exponents of a disordered system do not change, if 
the heat capacity of the pure system does not diverge. Later, a statement that the correlation 
length critical exponent of d-dimensional systems with such type of disorder should obey an 
inequality ly > 2/d was proven ( Chayes J.T., L. Chayes, Fisher D.S., T. Spenser, Phys. 
Rev. Lett., 1986, 57, 2999). For the extended, long-range correlated disorder the generalized 
Harris criterion holds (see the above mentioned paper of Weinrib and Halperin as well as D. 
Boyanovsky, Cardy J.L., Phys. Rev. B, 1983, 27, 6971). 

Y. Imry and S.-k Ma have shown (Phys. Rev. Lett., 1975, 35, 1399) that even arbitrary weak 
disorder coupled to the order parameter of continuous symmetry (i.e. for to ^ 2 vector model) 
destroys ferromagnetism at d < 4. For the random-field Ising model (to — 1) the lower critical 
dimension is d = 2. As I already have mentioned in the lectures, the isotropically distributed 
random axis destroys ferromagnetism at d < 4. One of the ways of showing this is to exploit 
the arguments similar to those of Imry and Ma for the random-field systems (Pelcovits R.A., 
Pytte E., Rudnick J., Phys. Rev. Lett., 1978, 40, 476). 

Last but not least, let me mention the Luck criterion for the connectivity disorder, as ex- 
plained in the lecture by Wolfhard Janke. 
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